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ABSTRACT 

We demonstrate that low resolution Call triplet (CaT) spectroscopic estimates of the 
overall metallicity ([Fe/H]) of individual Red Giant Branch (RGB) stars in two nearby 
dwarf spheroidal galaxies (dSphs) agree to ±0.1-0.2 dex with detailed high resolution 
spectroscopic determinations for the same stars over the range —2.5 < [Fe/H] < —0.5. 
For this study we used a sample of 129 stars observed in low and high resolution 
mode with VLT/FLAMES in the Sculptor and Fornax dSphs. We also present the 
data reduction steps we used in our low resolution analysis and show that the typical 
accuracy of our velocity and CaT [Fe/H] measurement is ^2 km s _1 and 0.1 dex 
respectively. We conclude that CaT- [Fe/H] relations calibrated on globular clusters 
can be applied with confidence to RGB stars in composite stellar populations over the 
range -2.5 < [Fe/H] < -0.5. 

Key words: techniques: spectroscopic - stars: abundances - galaxies: dwarf, star 
clusters 



1 INTRODUCTION 

An important aspect for a full understanding of galactic evo- 
lution is the metallicity distribution function of the stellar 
population with time. 

Carrying out detailed abundance analyses with high res- 
olution (HR) spectroscopy to trace the patterns that allow 
one to distinguish between the different galactic chemical 
enrichment processes is time consuming for large samples 
of individual stars in a galaxy. This is partly due to the 
observing time required, but also because of the complex 
data reduction and analysis necessary. Fortunately, there is 
an empirically developed, simply calibrated method avail- 
able which can make an efficient estimate of metallicity 
([Fc/H]) for individual Red Giant Branch (RGB) stars us- 
ing the strength of the Call triplet (CaT) lines at 8498, 
8542, 8662 A. This method was pioneer e d for use on indi- 
vidual stars by lArmandroff fc Da Costal l|l99lf) . It has the 
advantage that the lines are broad enough that they can be 
accurately meas ured with moderate spectral resolution (e.g, 
ICole et alJliuOi h 
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The CaT method is routinely used to estimate [Fe/H] 
for nearby resolved stellar systems and also provides an ac- 
curate radial velocity estimate. Both measurements are facil- 
itated by the strength of the CaT lines and by the generally 
red colours of the target stars. However, the CaT-derived 
abundances are empirically defined, with a poorly under- 
stood physical basis. Therefore it is important to check the 
results against HR spectroscopic (i.e. direct) measurements 
of [Fe/H] and other elements. The "classical" CaT calibra- 
tion is based on the use of globular cluster stars, all which 
are drawn from a single age and metallicity stellar popula- 
tion. The CaT equivalent widths are directly compared to 
HR spectroscopic measurements of [Fe/H] over a range of 
metallicity, and this comparison is used to define the rela- 
tion between CaT equivalent width (EW) and [Fe/H] for all 
observations taken with the same set up. This approach has 
been extensively tested for a large sample of globular clusters 
(Rutl edge et al. 1997a b, and references therein). However, 
globular clusters typically exhibit a constant [Ca/Fe] for a 
large range of [Fe/H]. This leads to uncertainty in the effect 
of varying [Ca/Fe] ratios such as is seen in the more complex 
stellar populations found in galaxies. Furthermore, stars in 
dwarf galaxies invariably cover a significant range of ages 
as well as metallicities. This mismatch in the properties of 
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calibrators and targets has led to suggestions that the CaT 
method may not be a very accurate indicator of [Fe/H] for 
more complex stellar populations, esp ecially in those c ases 
where [Ca/Fe] varies significantly fe.g. lPont et al.ll2004l ). 

In this paper we investigate the validity of the 
CaT method for complex stellar populations. We com- 
pare large samples of [Fe/H] measurements coming from 
VLT/FLAMES made using both the CaT method and di- 
rect HR spectroscopic measurements for the same stars in 
two nearby dwarf spheroidal galaxies (dSphs), Sculptor and 
Fornax, over a range of [Fe/H] and [Ca/H]. This is the first 
time such a detailed comparison has been made for stars out- 
side globular clusters. We also investigate the theoretically 
predicted behaviour of the CaT method for a range of stellar 
atmospheric parameters us i ng a grid of model atmosphere 
spectra from iMunari et all i|2005l ). 

The paper is organised as follows. In Section [2] we de- 
scribe the data reduction steps we use within the DART 
(Dwarf galaxy Abundances and Radial velocities Team) col- 
laboration to estimate EW and velocities from observations 
in the CaT region, as the accuracy with which this can be 
done clearly has important implications for the reliability of 
our conclusions for these galaxies. We also discuss the verifi- 
cation of the overall calibration and accuracy of the velocity 
and EW measurements by comparison of results from inde- 
pendent CaT observations and by comparing with theoreti- 
cal expectations based on signal-to-noise, resolution and line 
profile properties. In Section [3] we derive the standard CaT- 
[Fe/H] globular cluster calibration for low resolution (LR) 
VLT/FLAMES data. In Section Q] we compare the derived 
[Fe/H] from the CaT to the HR [Fe/H] for the Sculptor and 
Fornax dSphs. Finally, in Section [5] we discuss the uncer- 
tainties that come from using Call lines to derive an [Fe/H] 
abundance for stellar populations where the a-abundance 
varies and use a comparison with stellar model atmospheres 
to further investigate age, metallicity and a-abundance ef- 
fects. 



2 LOW RESOLUTION DATA REDUCTION & 
ANALYSIS 

The datasets presented here were collected between August 
2003 and November 2005. They consist of 15 pointings in 
the Sculptor dSph and 11 in the Fornax dSph spread over 
the galaxies (Fig. [1} . Some fields were observed with 1 hour 
exposure time, whilst other fields have repeated exposures of 
shorter integration time and two different plate set-ups, with 
the aim of testing the reliability of the derived velocities, 
EWs and the stability of the instrument. 

All the LR CaT observations were made using 
VLT/FLAMES in Medusa mode. This allows the simulta- 
neous allocation of up to 132 fibres, including dedicated 
sky fibres, over a 25' diameter field-of-view. We used the 
GIRAFFE low resolution grating (LR8), which covers the 
wavelength range from 8206A - 9400A, and gives a resolu- 
tion of R~6500. This allows the measurement of EW from 
CaT lines and also enables the derivation of velocities ac- 
curate to a Rifew km/s. This set-up was used as part of 
the DART programme to obtain spectra for several differ- 
ent fields in each of the Sculptor, Fornax and Sextans dSphs, 
and also for calibration purposes on a sample of 4 globular 
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Figure 1. Locations of the observed FLAMES targets 
at LR (squares) for t he Scu l ptor dSph (top), centred at 
01 ,1 00 m 09 s -33°42'30" llMated fi998h with position angle 99° 
lllrwin fc Ha tzidimitriou 1995), and for the Fornax dSph (bot- 
tom), centred at 02 fe 39 m 52 a -34°30'49" with position angle 
46.8° (Battaglia et al. 2006]). The ellip ses ind icate the tidal ra- 
dius ( Sculp tor: Irwin & Hatzidimitriou 1 19951 ; Fornax: Battaglia 
et al. 2006). The dashed circles show the FLAMES pointings. 
North is up and East is to the left. 



clusters: NGC104, NGC5904, NGC3201, NGC459(H, which 
cover the range -2.0 < [Fe/H] < - 0.7 on the CG97 scale 
(see Tabled]). 

The central fields in Sculptor and Fornax were also ob- 
served with a similar VLT/FLAMES set-up but at HR with 
R~20000, which facilitates direct measurement of individual 
lines, and hence direct abundance determination, of numer- 
ous elements (see Sect. [4]). 

Table[2]shows the journal of the VLT/FLAMES LR and 
HR observations we used for our analysis of the CaT- [Fe/H] 
calibration (Sects. [3] H}. 

The data were all initially reduced using the GIR- 



1 We acquired the data for NGC104 during an observing run in 
January 2005, whilst the data for the other globular clusters are 
from the ESO archive. 
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BLDRS^I pipeline pr ovided by the FL AMES consortium 
(Geneva Observatory. iBlecha et all 2003). This package pro- 
vided flat-fielding (including fringing removal), individual 
spectral extraction and accurate wavelength calibration, 
based on daytime calibration exposures. At the time we 
started this project no sky subtraction was available within 
this pipeline which led us to develop several further reduc- 
tion stages for the LR analysis. We describe the LR analysis 
in detail here. In Tables [3] and U we present the LR and HR 
results relevant for our analysis of the CaT-[Fe/H] calibra- 
tion (Sects. [3] H0). 

2.1 Sky subtraction and wavelength calibration 

An example of a spectrum produced by the GIRBLDRS 
pipeline is shown in Fig. The numerous skylines visible in 
this part of the spectrum not only serve as an independent 
check on the overall wavelength calibration, but also enable 
an update of the wavelength calibration of the individual 
spectra. 

Since the CaT lines only occupy a limited part of the 
LR spectra we optimised the sky subtraction and wavelength 
refinement for the range 8400A - 8750A . 

The first step is to combine all the sky spectra (typi- 
cally 10-20 sky fibres were allocated per field) using k-sigma 
clipping to remove spurious features and obtain an average 
sky spectrum. The result is then split into continuum and 
sky line components using an iterative k-sigma clipped non- 
linear filter (a combination of a median and a boxcar). The 
average "sky-line" spectrum is then used to define a sky-line 
template mask, in order to isolate those regions of the sky 
spectrum with significant features and mask out the remain- 
der. 

The processing of individual object spectra then pro- 
ceeds as follows: 

i. each spectrum is filtered as above to split the spectra 
into a line and continuum component, but this time addi- 
tionally masking out those regions affected by sky lines to 
allow a more accurate definition of the continuum; 

ii. the object line spectrum, which includes sky lines, 
is then cross-correlated with the masked line component of 
the average sky spectrum. This provides an accurate differ- 
ential wavelength update. The object line spectrum is then 
(rc)interpolated to be on the same wavelength scale as the 
average sky spectrum; 

iii. for the sky subtraction we compare the masked sky- 
line and object-line spectra, and find the optimum scale fac- 
tor and profile matching kernel that produce the minimum 
average absolute deviation (LI norm). This is applied to the 
line-only spectra. 

The optimum scaling factor, which in this case is cho- 
sen to minimise the LI norm rather than the commoner L2 
norm to reduce sensitivity to non-Gaussian outliers, is de- 
rived using a simple grid search with progressively finer step 
size. As noted previously, a mask is used to isolate the rel- 
evant regions of the sky spectrum to match to a template. 
By first removing the continua from both sky and object, 
more emphasis is placed on minimising the impact of sky 
line residuals. 

2 available at SourceForge, http://girbldrs.sourceforg.net/ 



iv. finally the object continuum is added back to the 
wavelength-updated object-line spectrum, the sky contin- 
uum is removed used the sky-line scaling factor and the 
sky-subtracted spectrum is saved for the next stage in the 
processing. Of course, implicit in the sky correction is the 
reasonable assumption that the derived scale factor for both 
sky lines and sky continuum is the same. 

The sky subtraction process described above involves 
some key components which are crucial to achieve good re- 
sults. Accurate wavelength registration is absolutely vital for 
good sky subtraction and is facilitated by the presence of co- 
pious numbers of strong sky lines. These sky lines are unre- 
solved at this resolution and at the signal-to-noise achieved 
(see later) which readily enables sub-km s _1 precision in 
wavelength alignment. This also has the added advantage of 
ensuring that systematic offsets due to wavelength calibra- 
tion for velocities from different observations are negligible. 

As a final step in this process, all the average sky spectra 
from each FLAMES observation, are cross-correlated with 
a chosen reference sky spectrum and used to put all the 
observations on the same internal system. This is done to 
avoid possible systematic differences between observations 
taken at different times. 

The effects of combining the sky spectra to form an av- 
erage sky and re-interpolating the object spectrum to this 
average wavelength system, almost invariably results in a 
slight mismatch between the spectral line profiles of the ob- 
ject and average sky spectra. This is circumvented by ap- 
plying a Hanning smoothing kernel to each in turn and 
finding which combination of smoothed and un-smoothed 
gives the best results (as determined by the optimum scale 
factor). More sophisti cated adaptive kernel matching (e.g., 
Alar d fc Luptonll 19981 ) is probably unwarranted in this case. 

2.2 Velocity and equivalent width estimation 

Our goal is to produce a robust automatic procedure that 
gives close to optimum results in terms of signal:to:noise and 
also produces minimal systematic bias as a function of EW. 

The first stage of the process is to estimate the con- 
tinuum in a similar manner to that described in the pre- 
vious section. Each spectrum is split into a line-only, and 
smoothly varying continuum-only component, using an iter- 
ative k-sigma clipped non-linear filter. This time the region 
around each CaT line is masked out, prior to filtering, to 
prevent the continuum tracking the wings of strong CaT 
lines. The effective scale length of the filtering is set to ~ 
15A which, in conjunction with the masking and iterative 
clipping, is sufficient to follow continuum trends without be- 
ing affected by the presence of strong lines. At the same time 
an estimate of the overall signal-to-noise in the continuum is 
made by measuring the median continuum level in the CaT 
region and the pixel-to-pixel noise covariance matrix in re- 
gions containing no lines. The latter is needed to correct 
the apparent random noise for the cumulative smoothing 
effects of spectral interpolation and resampling, which we 
find typically results in a factor of ~2 overall random noise 
reduction. 

After normalising by the computed continuum the ve- 
locity is estimated by cross-correlating each spectrum with a 
template. This template is constructed from a zero-level con- 
tinuum superimposed on three Gaussian absorption lines lo- 
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Figure 2. An example of a LR spectrum from FLAMES in the CaT wavelength range. In the upper panel is the unskysubtracted 
spectrum and in the lower panel is the result of automated sky subtraction. The CaT lines are marked. 



cated at the vacuum rest wavelengths of the CaT lines. Most 
of the weight in the least-squares fit comes from the core of 
the line which is sufficiently Gaussian-like in LR data, that a 
Gaussian fit provides an estimate with effectively a minimal 
rms error. Using more complex line profiles with correspond- 
ingly greater numbers of free parameters generally makes the 
rms error worse and is also more prone to wildly unstable so- 
lutions due to the inevitable presence of occasional artefacts 
in the data. The Gaussian line depths are scaled in the ratio 
3:5:4 to reflect the true relative strengths of CaT lines and 
all are set to have a full width at half maximum (FWHM) 
= 2.35A . An example of a continuum fitting to a Sculptor 
dSph K-giant spectrum together with the computed cross- 
correlation function and Gaussian fit to the peak is shown 
in Fig. H 

An accurate estimate of the position of the cross- 
correlation peak is made by fitting a Gaussian to a localised 
region around the peak. This velocity is then used to define 
the wavelength region around each CaT line to use for EW 
estimation. We estimate the EW in two ways. The first con- 
sists in simply summing the flux contained in a region cen- 
tred on each CaT line. After some trial and error we settled 
on a region 15A wide centred on each line as a reasonable 
trade-off between including all the line flux and minimising 
the noise. To derive the second estimate we fit individual 
unconstrained Gaussian functions to each CaT line over the 



same wavelength region (see Fig.|3J). This also allows a semi- 
independent check on the accuracy of the derived velocity by 
providing three separate velocity measures with associated 
errors. The weighted sum of these velocity errors provides 
the basic error estimate for the velocity derived from the 
cross-correlation. We prefer to use the latter method for the 
velocity estimate since it is effectively a constrained model 
fit. As a final step the derived velocity is corrected to the 
heliocentric system. 

The combined EWs for CaT lines #2 and #3 (A 854 2, 
A8662) for both the integral and Gaussian fits are then com- 
pared and used to compute an overall correction to the 
Gaussian fit. This is necessary since the real line profile is 
a complex function of many parameters, and in particular, 
the dampening wings visible in strong lines are distinctly 
non-Gaussian in appearance. This means that the observed 
CaT lines have non-Gaussian wings which are progressively 
more visible as the EW increases. To compensate for this we 
compare the ensemble Gaussian fit to the conventional EW 
integration for each dataset analysed, by finding the best fit 
slope that links the two sets. This correction is equivalent to 
computing the average overlap integral between the real line 
profiles and the Gaussian fits and is accomplished by mea- 
suring the gradient between the two as a function of EW. 
By using all the data to do this we introduce no additional 
rms error and remove the majority of the systematic bias as 
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Figure 3. Top: the upper spectrum shows the CaT region for a K-giant in the Sculptor dSph (V=17.7, S/N sa 40/A, derived [Fe/H]caT = 
— 1.50 dex), showing an example of automated continuum fitting (dashed line); the lower spectrum shows the residuals after continuum 
removal and subtracting Gaussian model fits to the three CaT lines. Bottom: the derived cross-correlation function and associated 
Gaussian fit around the peak region (± 50 km s" 1 about the peak). The fit is so good that it is indistinguishable from the observations. 



a function of EW. We also automatically take care of any 
generic spectrograph-induced line profile abnormalities. For 
the range of EWs considered here this method works reli- 
ably (see Fig. 2]). In deriving this correction we neglect line 
#1 (A849s) as it is the weakest and would only add noise to 
the determination. 

In the rest of our analysis we use the Gaussian-derived 
EW estimator since although the simple integral is unbiased 
it is also significantly noisier than a Gaussian fit. We show 
in the next section how the measurements of the continuum 
level and the random noise error in the continuum can be 
used to estimate errors in the EWs and velocities. 



2.3 Error bounds for velocities and equivalent 
widths 

For detailed abundance work most lines of interest are weak 
e.g. , EW » 100mA. Notable exceptions are the CaT lines 
which are generally heavily saturated and on the damped 
part of the curve-of-growth. Turbulence plus rotation of late- 
type giants typically only broaden the line profiles by « a 
few km s _1 , hence the profile of the CaT lines, which are 
completely unresolved at R=6500, are dominated by intrin- 
sic broadening due to saturation (typically FWHM 2A - 3A ) 



and to a lesser extent by the resolution of the spectrograph 
(FWHM wl.3A). 

Despite the dampening wings, to first order the LR CaT 
lines can be reasonably well approximated by Gaussian func- 
tions (e.g. see Fig. [3} and we can use this to gain some insight 
into the limiting factors that determine the accuracy of the 
velocity and EW measurements. 

Here we define the resolution as FWHM = 2.35cr where 
<t is the Gaussian profile equivalent scale parameter. Since 



where I p is the "peak" flux, this implies that line saturation 
(i.e. I p — C, where C is the continuum level/A ) occurs 
when EW w FWHM. The intrinsic FWHM of weak lines in 
these late-type giants is typically only a few km s _1 , i.e. 
lines with EWs of ~100mA and above are saturated. As 
noted previously, the CaT lines are heavily saturated and 
typically have EWs well above lA . 

To a reasonable approximation, the noise in the con- 
tinuum (a n per A ) due to sky plus object dominates, and 
over any individual line C can be taken to be constant over 
the lines of interest. Therefore, the EW and its error due to 
random noise, AEW, are given by 



EW 



H 

c 



AEW 



C 



S/N 
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Figure 4. A comparison of the EW obtained from profile fitting (EW fit) versus the EW from the direct integral of the lines (EW line), 
as an example of the relationship between the two 



where the total line flux is rj, S/N is the signal-to-noise per 
A and w is the effective width (A) the line is integrated over. 
For a fitted Gaussian profile 



AA 



FWHM 2 



4tt<7 « 1.5 FWHM 



hence 



AEW = 



VT5 FWHM 
S/N 



(2) 



(3) 



The FWHM of the line and continuum signal-to-noise are 
the primary abundance error drivers from a random noise 
point-of-view. For example, the two strongest CaT lines, 
■^8542 , A8662 , used in our analysis have FWHM at a resolution 
of R=6500 of 2A -3A , which for a continuum signal-to-noise 
of 10/ A implies a lower bound on the combined EW error 
of ^0.3A . 

In a similar way we can place constraints on the accu- 
racy of measuring velocities. For Gaussian-like line profiles, 
which are good approximations even for saturated lines like 
the CaT, the minimum variance bound on the error in the 
estimated line position A is given by 



~ a 2 alVW 
var\X\ = — 



(4) 



where the line flux is rj and where the noise in the contin- 
uum, a n (J A), dominates and can be take n to be a constant 
over the region of interest (see Irwin 1985, for more details). 
Rewriting this in terms of the FWHM of the line, the EW 
and the continuum signal-to-noise, leads to 



(5) 



EW S/N 

where all measurements are per A. 

As an example, for observations at R=6500, the FWHM 
of the strongest CaT line is typically between 2-3A , while 
the EW of these lines is typically ~2A ,implying accura- 
cies of wavelength centring at a continuum signal-to-noise 
of 10/Aof around 0.2Aper line, or equivalently ~5 km/s us- 
ing all three CaT lines. This is the minimum signal-to-noise 
we consider acceptable in our analysis and is achievable on 
VLT FLAMES for V=20 objects in 3600s of integration. 

2.4 Errors for velocity and equivalent width from 
repeated measurements 

The number of independent measurements for each dataset 
is 1740 for Sculptor and 1359 in Fornax. For Sculptor we 
have 464 stars observed only once, 428 stars observed twice, 
73 stars observed 3 times, 39 observed 4 times and 9 observed 
5 times. For Fornax we have 816 stars observed only once, 
209 stars with double, 23 stars with triple, 14 stars with 
quadruple measurements. 

We test for the reliability of our velocity and related 
errors by analysing the distribution of velocity differences 
from double measurements. The jth observed velocity Vji 
for a star i can be considered a random variable which fol- 
lows a Gaussian distribution centred around the true value 
«tfue,i with a dispersion given by the velocity error (jji. The 
difference between two repeated measurements vi t i and vi,%, 
Avi — vi,i — V2,i, is a random variable following a Gaussian 
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Figure 5. Comparison between velocity measurements for stars with double measurements in the Sculptor (left) and Fornax (right) 
dSphs. Panels a, b) Distribution of velocity differences for all the stars (dashed line, 428 stars for Sculptor and 209 for Fornax), and for 
the stars with S/N per A > 10 and estimated error in velocity < 5 km s _1 for each measurement (solid line, 203 stars for Sculptor and 
138 for Fornax). The weighte d mea n velocity, rms dispersion and scaled median absolute deviation from the median (1.48*MAD = a 
robust rms e.g. Hoaglin et al. ll983l) are: -0.5 ± 1.3 km s" 1 , 25.1 ± 1.0 km s _1 , 6.7 ± 0.2 km s" 1 (dashed line) and 0.1 ± 0.3 km s _1 , 
3.1 ± 0.4 km s" 1 , 3.1 ± 0.4 km s _1 (solid line) for Sculptor; 1.4 ± 2.0 km s" 1 , 27.7 ± 1.4 km s" 1 , 3.2 ± 0.2 km s" 1 (dashed line) and 
-0.01 ± 0.36 km s" 1 , 2.6 ± 0.4 km s" 1 , 2.4 ± 0.4 km s" 1 (solid line) for Fornax. Panels c, d) As above but now the velocity difference 
is normalised by the predicted error. With these S/N and velocity error cuts the measured error in the velocity distribution is very close 
to the expected unit variance Gaussian (standard deviation = 1.2 and MAD = 1.1 for Sculptor; standard deviation = 1.0 and MAD = 
0.9 for Fornax). Panels e, f) Comparison of velocities for stars with S/N per A > 10 and error in velocity < 5 km s _1 . 



distribution centred around zero and with dispersion given 
by <Ti = a\ j + ffji- Thus, if both velocities and their er- 
rors are correctly determined, the distribution of velocity 
differences Avi normalised by <7; should be a Gaussian with 
mean zero and dispersion unity. Figure [5] shows that if we 
take into consideration all the stars with repeated measure- 
ments both for Sculptor and Fornax the distribution has 



large tails of stars with Av/a > 3. If we restrict ourselves 
to the stars with S/N per A > 10 and error in velocity < 
5 km s _1 , then the resulting distribution is very close to 
a Gaussian. For Sculptor the resulting standard deviation 
is 1.2 and the scaled Median Absolute Deviation (MAD) is 
1.1, with just 5 stars with Av/a > 3 (2%). For Fornax the 
standard deviation is 1.0 and the scaled MAD is 0.9. 
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Figure 6. Comparison between summed EW (EW2+EW3) measurements of CaT lines from integrated flux (a,b) and Gaussian fitting 
(c, d, e, f) for stars with double measurements, S/N> 10 and error in velocity < 5 km s" 1 in the Sculptor (left, 203 stars) and Fornax 
(right, 138 stars) dSphs. The solid lines in panel (e, f) indicate the 1 and 3 <r region for an error in summed EW given by cr^w = 6/(5/ AQ 
(see text). Assuming this error, the weighted average, r.m.s. and MAD for the differences in T,W from Gaussian-derived estimator (panels 
c, d) is 0.04 ± 0.04 A , 0.30 ± 0.05 A , 0.32 ± 0.05 A for Sculptor and 0.07 ± 0.05 A , 0.36 ± 0.05 A , 0.34 ± 0.05 A for Fornax. 



Hereafter we consider a S/N per A > 10 and an error 
in velocity < 5 km s _1 as the minimum requirements for a 
reliable determination of velocity and EW. 

Figure |S^,,b,c,d shows the comparison of summed EW 
for the two strongest CaT lines, EW = EW2 + EW3, which 
we use as indicator of metallicity (see Sect.©, for the stars 
with double measurements in Sculptor and Fornax respec- 
tively. As mentioned in Sect. 12.21 the Gaussian-derived esti- 
mator is considerably less noisy than the integral estimator 
and thus we use this estimator in all our analysis. As an 
extra criteria to guarantee the reliability of the estimated 
EW measurements we also impose a cut-off in the difference 
between the Gaussian-derived and the integral estimators of 
less than 2 A . 

We derive the error in EW for single observations and 



how this varies with S/N from the comparison of EW for 
stars with double observations using Eq. ©. Figure [E£,f 
shows that, as expected, the comparison of EW improves 
for increasing S/N. We find that the random error in EW 
from repeated measurements is well represented by (Jew ~ 
6/(S/N), and hereafter we will adopt this formula as es- 
timate of the errors in the single measurements. This is a 
factor ~ 2 larger than we estimated from theoretical cal- 
culations; however a larger error is not surprising given the 
effect of all the steps involved in the data reduction (e.g. 
sky subtraction and continuum estimation). Such an error 
in EW results in a [Fe/H] error of ~ 0.1 dex at a S/N per 
A of 20 (see Sect.©. 

Finally, the measurements for stars with multiple obser- 
vations were combined, weighting them by their error, and 
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this results in 1013 distinct targets for Sculptor and 1063 for 
Fornax. The final sample was carefully checked to weed out 
any spurious objects (e.g., broken fibres, background galax- 
ies, foreground stars, etc). Excluding the spurious objects, 
and the objects that did not meet our S/N and velocity 
error criteria, our final sample of acceptable measurements 
consists of 648 stars in Sculptor and 944 in Fornax. Among 
these, 2 stars in Sculptor and 1 in Fornax did not meet our 
EW criteria. 

As an indication of the good quality of our data, the 
median S/N and median error in velocity are 32.1 and 1.6 
km s^ 1 , respectively, for the Sculptor dataset, and 24.3 and 
1.7 km s _1 for the Fornax dataset. 



3 THE STANDARD CAT CALIBRATION 
WITH VLT/FLAMES USING GLOBULAR 
CLUSTERS 

The next step is to transform the CaT EW into metal- 
licity, [Fe/H]. The dependence of CaT line strength on 
metallicity is theoretically difficult to understand, however 
it has been empirically proved by e xtensive calibration 
RGB stars in globular clusters llArmandroff fc Zinnl 

al ll99lf ). 
lation of 

CaT EW measure ments for individual RG B stars in glob- 
ular clusters, which iRutledge et~ai1 (|l997bh calibrated with 
HR metallicities, proving the CaT method to be reliable and 
accurate in the range —2 .1 < [Fe/H] < — 0.6. 

As summarised in Rutled ge et al. I (|l997bh the line 
strength index EW , which is a linear combination of the 
EW of individual CaT lines, depends on [Fe/H], the gravity 
logg and T c g. As \ogg and T e ff decrease going up the RGB, 
it is possible to remove the effect of the gravity and tem- 
perature taking into account the position of the star on the 
RG B with respect to the Horizonta l Branch (HB). 

lArmandroff fe Da Costal (|l99lft . using globular clusters, 
showed that this is most efficiently achieved when defining 
a "reduced equivalent width": 



In the literature there are several examples: IRutledge et al 
1 1997a! ) used a weighted sum of the 3 lines; iTolstov et al 



(2001 



using KAjfcS stars m globul ar clusters ( Arm andrott 
19881; lOlszewski et al.l[l99ll ; lArmandroff fc Da Coste 
Rutledge et al.l 1 1997af ) presented the largest compil; 



W' = EW + /3 X (V - Vhb) 



(6) 



where Vhb is the mean magnitude of the HB. The advantage 
of using V — Vhb over for instance the absolute I magnitude 
or the V — I colour is that the slope (3 is constant, since 
it does not vary with [Fe/H]. Using Vhb also removes any 
strong dependence on the distance and/or reddening. The 
"reduced equivalent width" W' is thus the CaT line strength 
index at the level of the HB. With the above definition, 
lArmandroff fe Da Costal (jlQQlh empirically proved that W' 
is directly correlated with [Fe /H], 

Using 52 globular clusters IRutledge et alj (|l997bh found 
P = 0.64±0.02Amag ~ 1 across the range —2.1 < [Fe/H] < — 
0.6. ICole et al.l (|2004r ) re-determined the value of f5 for their 
sample, which included globular and open clusters. They 
found (3 — 0.66 ± 0.03Amag _1 when usi ng only globu- 
lar clus ters, similar to the value found by IRutledge et al.l 
|l997br i, and (3 = 0.73 ± 0.04Amag _1 when including 
open clusters, which covered an higher metallicity range 
(-0.6 < [Fe/H] < -0.15). 

The first uncertainty in using the CaT calibration is of 
course which combination of the 3 CaT lines should be used. 



hereafter T01) excluded the first CaT line from the 



Cole et all (|2004l . hereafter C04) used an unweighted 



sum of the 3 lines. Such a choice usually depends on the 
quality of the data-set: in the case of limited signal-to-noise 
and with possible sky-line residual contamination, the weak- 
est line of the CaT is usually the least reliable and so the 
determination of its EW is often doubtful. 

These previous studies a ll calibrated the CaT W' on 
thelCarretta fc Grattonl |l997f ) scale (CG97), and lFriel et~ai] 
(12002) 1 for the open clusters, meaning that the relation be- 
tween [Fe/H] and W' has been derived assuming that the 
calibration clusters have the metallicities derived by CG97. 
To define their [Fe/H] scale, CG97 re-analysed high-quality 
EWs from different sources, using a homogeneous compila- 
tion of stellar atmosphere parameters, gf values and so on, 
making the largest and most self-consistent analysis of this 
kind. It is important to note that even the HR values of 
the [Fe/H] for globular clusters c an differ significantly from 
each other (see lPritzl et al.ll2005h . Obviously, choosing a dif- 
ferent [Fe/H] scale implies that the derived relations will be 
different. 

This variety of approaches and calibrations can lead to 
a degree of confusion when viewing the literature. It is not 
possible to find out a priori which way of summing CaT EWs 
and which calibration must be used; we thus test each of 
the mentioned approaches from the literature and see which 
one performs better for our globular clusters data-set. In 
addition, we derive our own calibration. 

Each of the three calibrations that we examined from 
the literature consists of three relations: 

1) a linear combination of the CaT lines EW: 

EW n97 = 0.5 EW X + EW 2 + 0.6 EWs (7) 

EWtoi = EW-2 + EWs (8) 

EW C0 4 = EW X + EW 2 + EW 3 (9) 

2) a relation for the "reduced equivalent width" : 
I¥ R97 = EVVW + 0.64(±0.02) x (V - Vhb ) (10) 
Wtoi = EIUtoi + 0.64(±0.02) x (V - Vhb) (11) 
Wcoi = £Wco4 + 0.73(±0.04) x (V — Vhb) (12) 

3) the calibration of the "reduced equivalent width" 
with [Fe/H]: 

[Fe/H]^ 7 97 = -2.66(±0.08) + 0.42(±0.02) x Wr 97 (13) 
[Fe/H]™ 7 = -2.66(±0.08) + 0.42(±0.02) x W^ 01 (14) 
[Fe/H]cG97 = -2.966(±0.032) + 0.362(±0.014) x Wc 04 (15) 



where R97 s t ands for IRutledge et~aT1 (Il997allbh . T01 for 
ITolstov et ail (|200ll ) and C04 for lCole et al l i|2004D . 

Figure [7] shows "EW versus (V— Vhb) for the 4 globular 
clusters, summing the CaT lines as in Eqs. (JT]), ([8}, (|9]). 
As the minimum S/N per A of these data is ~ 40 and the 
median is ~ 100 per A the errors in the summed EW are 
very small, <0.1 A . As consistency check we calculate the 
weighted average slope for each of the calibrations and find 
that they are consistent at the 1-cr level with the previous 
works (/3R97,thiswork = 0.59 ± 0.04 Amag -1 , /? T oi,thiswork = 
0.62 ± 0.03 Amag" 1 , ,3co4,thiswork = 0.79 ± 0.04 Amag" 1 ). 
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Figure 7. The CaT calibrations using globular clusters: EWs 
vs. (V — Vhb) f° r RGB stars in 4 globular clusters of different 
metallicities (asterisks: NGC104; squares: NGC5904; diamonds: 
NGC3201; triangles: NGC4590), using 3 different linear combi- 
nations of CaT lines (from top to bottom: R97, TOl, C04). The 
dotted lines show the relations in R97, TOl, C04, using the [Fe/H] 
published in CG97; the solid lines show our best-fitting relations, 
keeping the slope fixed within each calibration. The best-fitting 
metallicities are summarised in Table [T] 
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Figure 8. Relation between CG97 [Fe/H] and W' for the 4 cali- 
brating clusters using R97, TOl and C04 calibrations (from top to 
bottom). The solid lines show the best-fitting CG97 [Fe/H] and 
W' relation for each calibration. 



The metallicities derived from the best-fitting W' from 
an error- weighted linear fit of T.W versus (V — Vhb) are 
summarised in Table [T] In general the metallicities derived 
from these observations agree within la with the metallici- 
ties on the CG97 scale. The metallicities from the C04 cal- 
ibration seem to be systematically lower by ~ 0.1 dex. The 
TOl calibration appears to give the best performance. Thus, 
amongst the three relations in the literature, we will apply 
the TOl to our dSphs dataset. Figure [5] shows that the rela- 
tion between W' and [Fe/H]cG97 derived from the 4 globular 
clusters dataset is linear in each case. 

In order to derive our own calibration, we should repeat 
all the steps, i.e. derive a W' using the slope we find by 
fitting the summed EW versus (V — Vhb), and finding the 
best W'- [Fe/H] relation. However, as mentioned before, the 



Analysis and calibration of Call triplet spectroscopy 11 



slopes we find are consistent with those in the literature, and 
since the number of stars in most of our calibration globular 
clusters is not large, we prefer to use as relations those in 
the literature (Eqs. [T0lllllll2|) . We just repeat the last step 
and find the best-fitting W'-[Fe/H] relations by performing 
an error- weighted linear-fit. The relations we find are all 
consistent within 1-a with the relations in the literature. 
The calibration obtained using Eq. appears to give the 
best results and is the following: 

[ Fe /H]cG9 W 7° rk = (-2-81 ± 0.16) + (0.44 ± 0.04) x W4oi (16) 

In the following section we test to see if the [Fe/H]-W' 
relation derived for the globular cluster sample is reliable 
when applied to RGB field stars in galaxies. To do so, we 
apply Eqs. (|14[) and (|16[1 to the Sculptor and Fornax sam- 
ple and we compare the metallicities so derived to the HR 
metallicities. 



4 COMPARISON TO HIGH RESOLUTION 
METALLICITY MEASUREMENTS 

One major uncertainty in applying the CaT method to field 
stars in galaxies, for example in dSphs, is that the [Fe/H]- 
CaT W' relations have so far been calibrated exclusively on 
globular clusters. Some dSphs contain intermediate age and 
even young stellar populations and have a large spread in 
metallicity, whilst the above relations have been derived for 
single age stellar population older than 10 Gyr, over a rela- 
tively narrow metallicity range and also a very narrow range 
in [a/Fe], which is very different from composite stellar pop- 
ulations. Furthermore, in composite populations it is more 
difficult to assign a uni que magnitude for the HB, although 
ICole et all (|2000| . |2004| ) showed that the uncertainty due to 
this effect is ~ 0.05 dex, which is not significant compared 
to the intrinsic precision of the method (~ 0.1 dex). 

The only way to reliably test the standard globular clus- 
ter calibration for dSph field stars is to compare the [Fe/H] 
derived from CaT EWs to that obtained from direct mea- 
surements in HR observations of the same field stars. HR 
metallicities should be more accurate than CaT measure- 
ments because the iron abundance is not inferred from other 
elements but obtained by direct measurement of typically 
more than 60 separate Fe lines with two different ionisation 
states (i.e. Fel and Fell). 

Until now a comparison between HR and LR [Fe/H] 
has not been thoroughly made for dSph field stars, partly 
due to the lack of a large sample of overlapping measure- 
ments. Thanks to instruments like FLAMES, it is now pos- 
sible to get suitable comparison spectra for many objects at 
the same time. 

As part of the DART Large Program at ESO, HR 
FLAMES spectra have been taken for 93 probable Sculptor 
velocity member stars in the central regions of the Sculp- 
tor dSph (Hill et al. 2008, in preparation ) for which there 
is also LR CaT data (jTolstov et al.ll2004l ). A similar study 
was made of a central field in Fornax (Letarte et al. 2008, 
in preparati on; Letarte 2007), for which 36 stars overlap our 
LR sample (|Battaglia et al.ll2006l ). 

These observations consist of R~20000 resolution spec- 
tra of ~80 stars in the centre of both Fornax and Sculp- 
tor, covering ~60nm in three different FLAMES set-ups 



(534-562nm; 612-641nm; 631-670nm), and reaching typ- 
ical S/N of 30 per 0.05nm pixel. The chemical analy- 
sis of the sample was performed using OSMARCS one- 
dimensi onal stellar atmosphere model s in LTE (OSMARCS 
models, iGustafsson et all 120031 . 120071 ). an extension of the 
OSMARCS models referenced above (Pletz, private commu- 
nication), and a standard EW analysis. Stellar parameters 
were determined using a combination of photometric indices 
(V, I, J, H, K) and spectroscopic indicators (excitation and 
ionisation equilibrium). The results include the abundances 
of ~10 elements, including iron and calcium which are re- 
ported here for comparison to the LR results (see Table Q. 
Error bars on HR abundances indicated on the plots refer 
to the combined abundance measurement errors and propa- 
gated stellar parameters uncertainties. 

The detailed description of the data reduction and anal- 
ysis of the HR spectroscopic data can be found in a series of 
papers (Hill et al. 2008; Letarte et al. 2008; Letarte 2007). 
Also note th at the HR re sults have been put on to UVES 
system fe.g.. lLetartell2007l ). 

• CaT calibration using HR data: First we deter- 
mine the relation between [Fe/H]HR and CaT W' directly for 
the dSph data. We assumed Vhb = 20.13 for Sculptor and 
Vhb = 21.29 for Fornax, taken from llrwin fe Hatzidimitrioul 
(1995). Figure [9] shows the HR [Fe/H] of the overlapping 
stars between the HR and LR samples plotted against their 
reduced CaT equivalent width (Eg. 1 1 1 p . The best linear fit 
we obtain allowing for errors in both coordinates is 

[Fe/H] HR = -2.94(±0.04) + 0.49(±0.01)VF' (17) 

This calibration is consistent at the 1-a level with the cal- 
ibration derived in the previous section for 4 globular clus- 
ters, however there are some differences. For small W' the 
above relation predicts a [Fe/H] ~ 0.1 dex lower than from 
the globular cluster calibration given in Eq. (|16[1 , whilst the 
opposite happens for large W' . 

• Comparison between HR and LR results us- 
ing FLAMES globular cluster calibration: The tra- 
ditional globular cluster calibration of CaT is now applied 
to our LR CaT W' and the results compared to [Fe/HjHR- 
Figure \W\ shows the comparison between HR [Fe/H] and 
CaT [Fe/H] from our globular cluster calibration (Eq. I16[) . 
As also indicated in the previous figure, the two methods 
are generally in good agreement. The average difference is 
A [Fe/H] = [Fe/H] LR - [Fe/H] H R = -0.04 ± 0.02 dex and 
the spread is 0.17L0.02 dex, which is comparable with the 
measurement errors. We can thus apply the globular cluster 
calibration to dSph field stars with some confidence, between 
— 2.5< [Fe/H] < — 0.8 dex, where the relation between HR 
and LR data is linear. It is unclear if the comparison at 
[Fe/H]> —0.8 dex indicates a non-linearity in the relation- 
ship, as the appearance of "non-linearity" is given by ~ 4 
stars. To make a concrete statement requires more data at 
high [Fe/H]. 

The average difference calculated for the entire sample 
would suggest the absence of systematics in our evaluation of 
[Fe/H]LR. However when plotting A [Fe/H] versus [Fe/H]HR 
(bottom panel Fig. I10p a trend is visible, such that the LR 
[Fe/H] is overestimated of ~ 0.1 dex at [Fe/H] H R< - 2.2, 
and underestimated of ~0. 1-0.2 dex at the high [Fe/H] end, 
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Figure 9. Top: HR [Fe/H] versus CaT IV' for the RGB stars 
overlapping between the HR and LR sample for the Sculptor and 
Fornax dSphs (93 in Sculptor, squares, and 36 stars in Fornax, 
asterisks). The solid line represents the best error- weighted linear 
fit to the data (-2.94(±0.04) + 0.49(±0.01) W). Bottom: Resid- 
uals to the linear fit. The l<r y— errorbars are the errors in HR 
[Fe/H], The r.m.s. around the relation is 0.19 dex (MAD = 0.12 
dex), comparable to the measurement errors. 



at [Fc/HJhr > — 1.2, which instead suggests the presence of 
systematics in the [Fe/H]LR derivation. 

• Comparison between HR and LR results using 
globular cluster calibrations from the literature: As a 

comparison Fig. 1111 shows the CaT [Fe/H] de rived from the 
T01 c a libration (Eq. | 1 4 |). wh ich we used in Tolstoy et ail 
i|2004h . iBattaglia et all (|200rJ ) and lHelmi et all l|2006h . The 
trend with [Fe/H] is similar to the previous case, because of 
the similar slopes in the [Fe/H] versus W' relations. The 
systematic shifts are more enhanced with respect to the 
globular cluster calibration at the low [Fe/H] end, whilst 
the performance at the high [Fe/H] end is better than for 
Eq. (0. 



Both the relations in the literature and our FLAMES globu- 
lar clusters calibration are in agreement with the HR data in 
the range -2.5 < [Fe/H] < -0.8 dex. Given the slight differ- 
ences between the relations, the choice of one over the other 
does not affect considerably the derived metallicities. In the 
remainder of the paper we follow the standard approach 
and use our own globular cluster calibration (Eq. I16p . We 
also checked that [Fe/H] derived from our CaT calibration 
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Figure 10. Top: Comparison between HR and CaT [Fe/H] for 
the Sculptor (squares) and Fornax (asterisks) dSphs. The solid 
line indicates the one-to-one relation. The CaT metallicities are 
derived using Eq. 116> . Bottom: The above comparison plotted 
as a difference. The diamonds with errorbars show the weighted 
average (and dispersion) of [Fe/H]LR — [Fe/H]n R , in bins of 0.3 
dex. The average difference is [Fe/H] LR -[Fe/H] H R = -0.04±0.02 
dex. The r.m.s. scatter around the one-to-one relation is 0.17 dex 
(MAD = 0.12 dex). For the stars with [Fe/H] HR > -0.8 the 
scatter and MAD are 0.2 dex and 0.21 dex, respectively. 



is in a greement with other derivations, such as in IPont et al.l 
d2004h . 



5 POSSIBLE SOURCES OF UNCERTAINTY IN 
THE CAT METALLICITY 
DETERMINATION 

Although [Fe/H] derived from the CaT method is in good 
agreement with the derivation from HR measurements, it 
would be interesting to understand the causes of the appar- 
ent systematics. 

As mentioned in Sect. [4] the traditional CaT calibration 
has a number of implicit assumptions in it: 1) [Ca/H] does 
not play a role in determining W , and thus one can exclude 
it from the calibration; 2) alternatively, if it plays a role, 
one assumes that globular clusters and dSphs have similar 
[Ca/Fe]; 3) the CaT EW is a better estimator of [Fe/H] than 
[Ca/H]; 4) the effect of age differences between dSphs and 
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Figure 11. Top: As previous figure, but the CaT metallicities are 
derived from Eq. H14I I. Though the average difference is very small, 
[Fe/H] LR - [Fe/H] H r =0.04±0.02 dex, the bottom panel shows 
that there is a gradient in [Fo/HJlr — [Fc/H]hr with [Fe/H], more 
enhanced than in the previous case. The r.m.s. scatter around the 
one-to-one relation is 0.17 dex (MAD = 0.12 dex). 



globular clusters stars can be neglected. In the following we 
explore the validity of these assumptions. 



5.1 CaT EW as [Ca/H] estimator 

5.1.1 Calibration from globular clusters 

To test assumption 3) we derive a relation between W' and 
[Ca/H] using the globular clusters as calibrators, and we 
apply it to our sample of dSphs. If W' traces [Ca/H], then 
we should expect a good one-to-one relation between the 
Ca abundance derived from the globular clusters calibration 
and the Ca abundance derived from the direct HR measure- 
ments. 

In order to derive the [Ca/H] values needed for the 
globular cluster W'-[Ca/H] relation we assumed an aver- 
age value of [Ca/Fe]=0.235 dex for the 4 globular clus- 
ters (see Table [T] for the individual values), and we used 
[Ca/H] = [Ca/Fe] + [Fe/H], where [Fe/H] are the individual 
values from CG97 (see Table [TJ. 

The best-fitting relation we obtain from the globular 
clusters is 




[Ca/H] c 



-2.57(±0.18) +0.44(±0.05)dex/A W' . (18) 
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Figure 12. [Ca/H]c aT from Eq. (TTHJ versus [Ca/H] H R for the 
Sculptor (squares) and Fornax (asterisks) dSphs. The solid line 
shows the one-to-one relation. 



Using the individual [Ca/Fe] values (instead of the average 
[Ca/Fe] =0.235 dex) does not change the best- fitting rela- 
tion, it just increases the resulting minimum \ 2 value of 
the fit. Figure [12] shows the comparison between [Ca/H]caT 
and the corresponding [Ca/H] derived from HR measure- 
ments for the Sculptor and Fornax dSphfl The agreement 
between the two is not as good as for [Fe/H]: the scatter is 
much larger and [Ca/H]caT is in general overestimated, es- 
pecially at [Ca/H] > — 1.3. The comparison gets worse if we 
approximate the trend of [Ca/Fe] with [Fe/H] with a func- 
tion, to circumvent the effect of the large errors in [Ca/Fe]. 

Thus the Ca abundance cannot be the only factor that 
drives the W' . 



5.1.2 Calibration from the literature 

We compared our results to those of iBosler et all (|2006T) . 
They used 15 Galactic globular clusters and 6 open clusters, 
covering the range 2 Gyr < age < 14 Gyr and —2.2 < 
[Ca/H] < +0.2, to derive a calibration between CaT W' 
and Ca abundance, the idea being that using W 1 as tracer 
of [Ca/H] instead of [Fe/H] one can apply the calibration to 
systems such as globular cluster or dSphs, without worrying 
about the different [Ca/Fe] trends. The relation they derive 
is: 

[Ca/H] B osi cr = -2.778(±0.061)+0.470(±0.016)dex/A W'{1Q) 

and they apply it to the Leo I and Leo II dSphs. This relation 
is consistent with our Eq. (|18[) . an d this results in a poor 
comparison between [Ca/H]caT and [Ca/H]HR- 

5.2 Exploring the effect of [Ca/Fe] 

Another of the uncertainties in applying the CaT method to 
composite stellar populations (i.e. galaxies) is the varied and 

3 Out of the 93 stars in the overlapping sample in Sculptor, 4 
had too low signal-to-noise to allow a determination of the HR 
[Ca/H] 
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Figure 13. HR [Ca/Fe] versus HR [Fc/H] for the Sculptor 
(squares) and Fornax (asterisks) dSphs. The triangles show the 
values for the calibration globular clusters (Table [TJ- The dia- 
monds with errorbars show a weighted average and dispersion of 
[Ca/Fe] in [Fe/H] bins. 




Figure 14. Observed CaT W' for our globular cluster sample ver- 
sus CaT W derived from a linear combination of Fe and Ca abun- 
dances (Eq. |20> . calibrated on the dSphs sample (squares) and 
from the linear dependence on Fe abundance given by Eq. ( jlTI 
(asterisks). The solid line indicates the one-to-one relation. 



extended star formation histories of these systems, which re- 
sults in a range of [Ca/Fe] values and a trend with [Fe/H]. 
Figure [13] shows that our calibrating globular clusters have 
an almost constant [Ca/Fe] trend with [Fe/H] (except for 
NGC3201), whilst both for the Sculptor and Fornax dSph 
[Ca/Fe] decreases with [Fe/H]. At [Fe/H] < - 1.1 there is 
overlap between the [Ca/Fe] values for the globular clusters 
and the dSphs, although [Ca/Fe] in the globular clusters is 
larger than the average for the dSphs at the same [Fe/H]; 
at [Fe/H]~ —1.1 the [Ca/Fe] value for the globular clusters 
is ~ 0.2 dex larger than the average value for Sculptor and 
at [Fe/H]= —0.7 is ~ 0.5 larger than the average value for 
the Fornax dSph. This results in a smaller [Ca/H] abun- 
dance for globular clusters with respect to our dSphs sam- 
ple at the low [Fe/H] end, and we have the opposite effect 
at the high [Fe/H] end. If the W' is an increasing function 
of both [Fe/H] and [Ca/H] abundances, to neglect the effect 
of [Ca/H] in our globular cluster calibration would result in 
an overestimated [Fe/H]LR with respect to [Fe/H]HR in the 
region where [Ca/H]ds P h > [Ca/H]Gc; instead we would un- 
derestimate [Fe/H]LR with respect to [Fe/H]HR in the region 
where [Ca/H]ds P h < [Ca/H]GC- This goes in the same direc- 
tion to what we see in Fig. 1101 This suggests that [Ca/H] 
might also play a role in determining the CaT W' . 

It is not obvious how to quantify this effect as it is 
not known how the CaT W' depends on [Ca/H] (linearly, 
quadratically, etc). We tested for a simple dependence of 
this kind: 

W' = a[Fe/H] + 6[Ca/H] + c (20) 

where a,b and c are constants. We derive a = 1.58 ± 0.09, 
b = 0.36 ± 0.13, c = 5.85 ± 0.07 by fitting the observed 
values of W' from the dSphs dataset as a function of the 
corresponding [Fe/H] and [Ca/H] from HR measurements. 
This relation, which suggests that [Ca/H] does have an im- 
pact on W , but it is less significant than [Fe/H], agrees well 
with the observed values of W' when applied to the globu- 
lar cluster data (Fig. I14[) . although it does not remove the 
trend altogether. The importance of [Ca/H] in driving W' is 
anyway unclear since the comparison between the observed 
W' for the globular clusters and the W' predicted by ap- 



plying Eq. (|17p - therefore neglecting the effect of [Ca/H] - 
improves slightly at the low W' end while gets slightly worse 
at intermediate W' . 

This analysis is clearly not exhaustive, but indicative 
of the effects of the obvious sources of uncertainty (see also 
ICenarro et al. 2002). Since we are interested in understand- 
ing the effect on the [Fe/H] determination, it is important to 
notice that notwithstanding the large difference in [Ca/Fe] 
between calibrating globular clusters and dSphs at the high 
[Fe/H] end (~ 0.5 dex at [Fe/H]- 0.7) and the large dif- 
ference in age (Fornax stars at ~ —0.7 have an age of 3-6 
Gyr) the error in estimating [Fe/H] using globular clusters 
as calibrators is just 0.1-0.2 dex. 

We point out that in this respect the Fornax dSph is 
one of the most extreme cases as the large majority of dSphs 
do not show such an extended star formation history, and 
so the difference in [Ca/Fe] (and in age) with respect to the 
calibrating globular clusters will typically not be as extreme. 

5.3 On age effects 

As shown in iPont et al.l <|2004l ), for metallicities < —0.5 dex, 
a relative age difference of ~12 Gyr with respect to calibrat- 
ing globular clusters results in an underestimate in metal- 
licity of ~0.4 dex (their fig. 16, middle panel). This error 
decreases substantially for smaller age differences. For in- 
stance, for a metallicity ~ —0.5, the error induced by us- 
ing a calibrating cluster 12.6 Gyr old for a star which is 4 
Gyr old is ~ ± 0.2 dex. At lower metallicities this effect is 
small er. The stars in the Scul ptor dSph have ages > 8 Gyr 
(e.g., iMonkiewicz et all [i~999) and [Fe/H] < - 1, thus the 
age effect will be negligible. The bulk of stars in the For- 
nax dSph have an age betw een 3-6 Gyr old (e.g.. |Pont et al.l 
2004; Battagli a et al.1 [20061 ); only a small fraction of stars 
have younger agefl, thus the age effect could result in an 
underestimate of ~0.1-0.15 dex from the CaT measurement 
for 3-6 Gyr old stars. As mentioned above, most of the dSphs 

4 these spectroscopic samples of RGB stars only contain stars 
older than 1-2 Gyr 
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do not show star formation histories as extended as Fornax, 
thus the age effect should typically be even smaller. 



5.4 Synthetic stellar atmospheres 

Another way of approaching these issues is to model the ex- 
pected behaviour of the CaT lines using stellar atmospheres 
and synthetic spectra. It is perhaps not surprising that these 
lines exhibit a strong dependence on global metallicity, best 
represented by [Fe/H], rather than on the abundance of cal- 
cium itself, given the strength of the features. With EWs 
ranging from ~500mA to several A , the CaT lines are not 
only saturated but are truly in the strong line regime of line 
formation, when a significant part of the flux is found in 
the pressure-broadened wings rather than the core of the 
line. This effectively means that electronic pressure (P e ) 
is driving the strength of the line rather than the usual 
combination of temperature and elemental abundance that 
shape weak lines. This can explain a fundamental charac- 
teristic of the CaT: since many metals contribute to P e , the 
CaT lines become sensitive to the global metallicity, [Fe/H], 
rather than calcium abundance alone, through the pressure- 
dependent wings. In fact, natural broadening dominates the 
wings of CaT lines, so that line strength will increase with 
decreasing electronic pressure (roughly as 1/P e ). This ex- 
plains why the CaT line strengths increase with increasing 
luminosity (i.e. decreasing gravity and hence pressure). The 
metallicity dependence of P e will therefore also contribute to 
shaping the CaT lines. Finally, global metallicity also plays 
a role in changing the blanketing properties by numerous 
small metallic absorption lines, both in the wavelength re- 
gions used to define the continuum and in the CaT wings 
themselves. It is therefore difficult to speculate theoretically 
how the CaT W' should behave upon varying each stellar pa- 
rameter and we have therefore used synthetic stellar spectra, 
"observed" with an approach mimicking our observational 
procedure. 

We have checked the validity of calibrating [Fe/H] to 
the width of the CaT lines by using synthetic spectra of 
RGB stars. |m unari et alj (|2005l ) published a large grid of 
synthetic spectra covering the CaT region at various res- 
olutions. Since our FLAMES LR resolution of R « 6500 
did not correspond to the available resolutions we rebinned 
the Munari R = 20000 resolution spectra to the FLAMES 
LR resolution. The model atmosphere spectra cover a range 
of stellar parameters including logg, T e g, [M/H] and also 
include models computed with [a/Fe]=0.0 and +0.4. The 
CaT EWs of a representative sample of model atmosphere 
spectra covering the range of stellar RGB parameters en- 
countered in dSphs were measured in a similar manner to 
those of our observed LR spectra. 

In particular, we measured the EW of a set of syn- 
thetic stars taken along the upper RG B of a series of 12Gy r 
isochrones with -2.5 < [Fe/H] < +0.0 (|Bertelli et al.lll994h . 
We then derived the corresponding [Fe/H] using our CaT 
calibration based on globular clusters (Eq. I16|l . The V mag- 
nitude for each synthetic star was read off the isochrone, 
while the Vhb was assumed [Fe/H] dependent, o f the form 
Vhb = 1.17 + 0.39X [Fe/H] (|Ortolani et al.lll995l , and refer- 
ences therein, which also fit well the same isochrones). The 
resulting (V — Vhb) range from 0.5 to 3.5 (RGB tip), and 



are very similar to the range of luminosities of our targets 
in Sculptor. 

The CaT EWs were measured in both a-poor and a- 
rich models, and Fig. H5l shows how remarkably well the re- 
sulting metallicities from the CaT calibration ([Fe/H] CaT) 
agree with the input model metallicity ([Fe/H]model). Our 
calibration on the a-rich synthetic spectra overestimates the 
[Fe/H]model by only +0.19 dex with a negligible disper- 
sion and no particular trend, while our [Fe/H]CaT of the 
a-poor spectra underestimates [Fe/H]model by —0.38 dex 
in the mean. This is expected, as our calibration is based on 
globular clusters, which have in the mean [ q/Fe]~ +0.25, 
and a re therefore closer to the a-rich grid of iMunari et al.l 
(2005). Moreover, the relation between the CaT W' and 
[Fe/H] model is linear, with a slope (for a-rich) matching 
exactly the slope of our calibration. This reassuring result 
supports the soundness of using simple CaT (linear) calibra- 
tions to derive metallicities. 

We also examined the effect of using the same calibra- 
tion for synthetic stars along younger RGB isochrones of 
2.5Gyrs (similar to the mean age of the stars in our For- 
nax sample). The main effect is the luminosity increase: 
at a given metallicity, the RGB isochrones of intermediate 
and old ages overlap almost perfectly in temperature and 
gravity, but younger RGB stars appear more luminous than 
old stars of the same temperature and gravity (because a 
younger star with the same gravity and temperature will 
be of higher mass, and hence more luminous). Between a 
12 and 2.5 Gyrs isochrone, the luminosity at a given gravity 
increases by ~0.5 mag, translating into a — 0.2dex metal- 
licity decrease once the CaT calibration is applied, as illus- 
trated in Fig. 1151 This age effect is small, and comparable 
to the shallow slope observed in Fig. \W\ (lower panel) where 
our younger metal-rich stars (Fornax) tend to give lower 
[Fe/H]CaT abundances than expected. 

As a final comment on this comparison, we would like 
to stress that the CaT behaviour as a function of luminos- 
ity and metallicity assumed in this work is fully consistent 
with what is expected from synthetic spectra. The synthetic 
spectra are however probably not precise enough to further 
constrain the CaT behaviour, as the core of these very strong 
lines is almost certainly not well modelled by the stellar at- 
mosphere models, as the cores of the lines form close to the 
stellar surface, a region which is problematic to model, and 
may even include a possible contribution from the stellar 
chromosphere. 



6 SUMMARY AND CONCLUSIONS 

We described the data reduction steps we use within the 
DART collaboration to estimate velocities and CaT metal- 
licities for LR data for RGB stars in dSphs. We showed that 
we obtain accurate velocities and [Fe/H] measurements, with 
internal error in velocity ~ 2 km s -1 and in [Fe/H] ~ 0.1 dex 
at a S/N per A of 20. 

We used 4 Galactic globular clusters observed with 
VLT/FLAMES in the CaT region to test the performance 
of several CaT W'-[Fe/H] relations existing in the liter- 
ature. We also derived the best calibration from these 
globular cluster data. The relation here derived is con- 
sistent at the 1-cr level with the calibration derived in 
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Figure 15. Left: Relation between the W' and the input metallicity of the synthetic spectra of the lMunari et aD l|2005h grid ([Fe/H]model) 
along a 12Gyr isochrone. Filled squares are the models with [a/Fe]=0.0 and filled triangles are the models with [a/Fe]=0.4. The thick 
black line shows our calibration. The open symbols show the effect of changing the age of the isochrone from 12 to 2.5Gyrs for metallicitics 
[Fe/H]> — 1. Right: Comparison between the metallicity obtained using our CaT calibration on the synthetic spectra ([Fe/H]CaT) and 
the input metallicity of the model ([Fe/H]model). Dotted lines show the mean offsets of the [a/Fe]=0.0 models and [«/Fe]=0.4 models 
with respect to the one-to-one relation (solid line). 



Tolstoy et all (l200lh. which we used inlTolstov et all l|2004h . 
Battaglia et al.l l|2006l ) and lHelmi et alj (|2006l ~ 



We used a sample of 93 and 36 RGB stars in the Sculp- 
tor and Fornax dSphs, respectively, overlapping between our 
LR and HR VLT/FLAMES observations, to test the globu- 
lar clusters CaT calibration. This is the first time that the 
CaT calibration is tested on field stars in galaxies. We find a 
good agreement between the metallicities derived with these 
two methods. However, a systematic trend is present with 
[Fe/H], such that using the globular cluster calibration de- 
rived in this work the [Fe/H] measurement from CaT is 
overestimated of ~0.1 dex at [Fe/H]HR < —2.2, whilst at 
[Fe /H] hr > -1.2 it is underestimated of ~0.1-0.2 dex. No 
clear systematic trend is instead derived from our data for 
[Fe/H]HR > —0.8. In order to understand this systematic ef- 
fect, we explored the possible contribution of Ca abundance 
to the calibration, and showed that there are indications 
that it might well affect the CaT W', although much less 
than [Fe/H]. From our dataset we also show that, contrary 
to previous claims, it is not advisable to use the CaT W' as 
a linear indicator of [Ca/H]. 

Finally we investigated the effect of varying stellar at- 
mosphere parameters on the CaT method b y analysing a 
large sample of model atmosphere spectra (|Munari et al.l 
2005). We a gain demonstrated that the CaT method is (sur- 
prisingly) robust to the usual combination of age, metallic- 
ities and [a/Fe] variations seen in nearby dSphs. 

From our analysis we see that even for large differences 
in [Ca/Fe] between calibrating globular clusters and our 
sample of dSphs (~ 0.5 dex at [Fe/H]~ —0.7) and large 
difference in ages (at [Fe/H] ~ —0.7 Fornax stars are ~ 10 
Gyr younger than globular clusters stars) the error in esti- 
mating [Fe/H] using globular clusters as calibrators is just 
0.1-0.2 dex. The Fornax dSph is likely to represent the most 



extreme case as it has had one of the most extended star 
formation histories among the dSphs in the Local Group. 

We conclude that CaT- [Fe/H] relations calibrated on 
globular clusters can be applied with confidence to RGB 
stars in composite stellar populations such as galaxies, at 
least in the [Fe/H] range probed by the above analyses, 
-2.5 < [Fe/H] < -0.5. Hence, the CaT method provides 
a good indicator of the overall metallicity of resolved stars. 

This has implications for the efficiency with which we 
can obtain metallicity distribution functions of nearby re- 
solved galaxies. The HR data collected in this paper required 
more than 6 nights of VLT observing time for ~ 150 spectra, 
whereas ~ 120 CaT spectra were obtained in one hour VLT 
observing time. 
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Cluster 


RA(J2000) 


DEC(J2000) 


V HB 


v sys [km s 1 ] 


[Fe/H]cG97 


[Fe/H] Rg7 


[Fc/H]toi 


[Fc/H]co4 


[Ca/Fe] 


NGC104 


00 h 24 m 05.2 s 


-72°04'51" 


14.06 


-14.4±0.7 


-0.70±0.03 


-0.78±0.05 


-0.68±0.06 


-0.82±0.06 


0.30 


NGC5904 


15 h 18 m 33.8 s 


+02°04'58" 


15.07 


53.2±0.7 


-1.11±0.03 


-1.13±Q.07 


-1.02±0.08 


-1.19±0.06 


0.21 


NGC3201 


10 h 17 m 36.8 s 


-46° 24' 40" 


14.77 


497.5±1.2 


-1.23±0.05 


-1.27±0.07 


-1.17±0.06 


-1.36±0.07 


0.11 


NGC4590 


12 fc 39 m 28.0 s 


-26°44'34" 


15.68 


-98.3±1.7 


-1.99±0.06 


-1.94±0.06 


-1.90±0.07 


-2.11±0.06 


0.32 



Table 1. The globular clusters used in this work. Positi ons and Vhb are fro m iHarrij dl996t ). The systemic velocity v sys was derived 
in this work. [Fe/H]cG97 i s the value of metallicity from [Carretta j z, Grattonl fll997f). [Fe /H]jj97, [Fe/H l xoi, [ F e/Hlnp4 are the v alues 
of the metallicity from this work, using the calibrations from iRutledge et al.l dl997al lbl), iTolstov et al.l ll2001h . ICole et al] d2004T) . i.e. 
Eqs. <13t . 114t . 115t . The errors in the derived metallicities were obtained from the stan dard deviation of the values of reduced equivalent 
width for each gl obular cluster. The [Ca/Fe] values are from ICarretta et al. | <2004h for 47 Tuc and all the other ones are listed in 
iKoch et al.l ||2006| '). 



Object 


RA(J2000) 


DEC(J2000) 


Date of observation [YY-MM-DD] 


cxptime [s] 


grating 


Sculptor 


01 


00 


03 


-33 


11 


30 


2003-08-20 


4200 


HR13 


Sculptor 


01 


00 


03 


-33 


41 


30 


2003-08-21 


3x3600, 4700 


HR14 


Sculptor 


01 


00 


03 


-33 


11 


30 


2003-08-22 


2x3600 


HR13 


Sculptor 


01 


00 


03 


-33 


11 


30 


2003-08-22 


4500 


HR14 


Sculptor 


01 


00 


03 


-33 


41 


30 


2003-08-23 


5400 


HR14 


Sculptor 


01 


00 


03 


-33 


41 


30 


2003-08-24 


2x3600 


HR10 


Sculptor 


01 


00 


03 


-33 


11 


30 


2003-08-25 


5400 


HR10 


Sculptor 


01 


00 


03 


-33 


11 


30 


2003-08-28 


3600 


HR10 


Sculptor 


00 


58 


22 


-33 


42 


36 


2003-09-29 


2x1800 


LR8 


Sculptor 


01 


01 


44 


-33 


43 


01 


2003-10-01 


2x1800 


LR8 


Sculptor 


01 


00 


03 


-33 


11 


29 


2004-09-10 


3600 


LR8 


Sculptor 


01 


01 


11 


-33 


56 


50 


2005-11-07 


3600 


LR8 


Sculptor 


00 


58 


59 


-33 


27 


04 


2005-11-08 


2x2700 


LR8 


Fornax 


02 


39 


30 


-34 


28 


39 


2003-09-29 


2x3600 


HR13 


Fornax 


02 


39 


30 


-34 


28 


39 


2003-09-29 


2x3100 


HR14 


Fornax 


02 


39 


30 


-34 


28 


39 


2003-09-29 


2x3600 


HR13 


Fornax 


02 


39 


30 


-34 


28 


39 


2003-09-30 


3300, 2x5400 


HR14 


Fornax 


02 


39 


30 


-34 


28 


39 


2003-09-30 


3600 


HR10 


Fornax 


02 


39 


30 


-34 


28 


39 


2003-10-01 


3x3600 


HR10 


Fornax 


02 


39 


30 


-34 


28 


39 


2003-10-01 


3300, 3600 


HR14 


Fornax 


02 


39 


30 


-34 


28 


39 


2004-01-14 


3600 


HR10 


Fornax 


02 


39 


30 


-34 


28 


39 


2004-01-15 


3600 


HR13 


Fornax 


02 


39 


30 


-34 


28 


39 


2004-01-19 


2x3600 


HR13 


Fornax 


02 


39 


30 


-34 


28 


39 


2004-01-20 


3600 


HR13 


Fornax 


02 


39 


30 


-34 


28 


39 


2004-01-21 


2x3600 


HR14 


Fornax 


02 


39 


30 


-34 


28 


39 


2004-01-22 


3600,3900 


HR14 


Fornax 


02 


39 


30 


-34 


28 


39 


2004-01-23 


3600 


HR10 


Fornax 


02 


39 


30 


-34 


28 


39 


2004-01-23 


3600 


HR14 


Fornax 


02 


39 


30 


-34 


28 


39 


2004-01-24 


3600 


HR10 


Fornax 


02 


39 


30 


-34 


28 


39 


2004-01-26 


3600 


HR10 


Fornax 


02 


38 


56 


-34 


24 


48 


2005-11-07 


3600 


LR8 


Fornax 


02 


39 


36 


-34 


40 


51 


2003-08-24 


1500,1800 


LR8 


Fornax 


02 


39 


58 


-34 


15 


28 


2003-08-26 


2400 


LR8 


NGC104 


00 


22 


15 


-72 


04 


45 


2005-01-11 


3x600 


LR8 


NGC3201 


10 


17 


36 


-46 


24 


40 


2003-02-22 


300, 327 


LR8 


NGC3201 


10 


17 


36 


-46 


24 


40 


2003-03-06 


2x600 


LR8 


NGC4590 


12 


39 


28 


-26 


14 


33 


2003-03-06 


3x600 


LR8 


NGC5904 


15 


18 


03 


02 06 40 


2003-03-05 


6x600 


LR8 



Table 2. Table of VLT/FLAMES observations for the Sculptor and Fornax dSphs and the globular clusters NGC104, NGC3201, 
NGC4590, NGC5904. We used these observations for the analysis of the CaT-[Fe/H] calibration. 
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Object Star ID 


RA(J2000) DEC(J2000) 


V 


eV 




cv hcl SWtoi 


eXWroi [Fc/H]lr, c[Fc/H] lr 


NGC104 GB07-TUC1002 





23 


41.71 


-72 5 2.6 


13.052 


0.002 


-8.55 


0.38 


5.209 


0.055 


-0.80 


0.02 


NGC104 GB07-TUC1510 





24 


25.01 


-72 2 45.7 


12.652 


0.001 


-6.93 


0.57 


5.729 


0.039 


-0.69 


0.02 


NGC3201 GB07-N3201-S60 


10 


17 


42.84 


-46 22 48.5 


13.601 


0.002 


495.55 


0.65 


4.144 


0.055 


-1.32 


0.02 


NGC3201 GB07-N3201-S74 


10 


17 


48.53 


-46 23 39.3 


14.602 


0.002 


498.98 


0.63 


3.636 


0.061 


-1.26 


0.03 


NGC4590 GB07-M68-S99 


12 


39 


40.21 


-26 47 37.5 


15.186 


0.004 


-100.62 


1.80 


2.021 


0.132 


-2.06 


0.06 


NGC4590 GB07-M68-S144 


12 


39 


32.06 


-26 47 53.6 


15.291 


0.004 


-99.36 


1.92 


2.139 


0.139 


-1.98 


0.06 


NGC5904 GB07-M5-S23 


15 


17 


55.26 


2 7 25.6 


14.544 


0.002 


55.11 


1.08 


4.175 


0.036 


-1.12 


0.02 


NGC5904 GB07-M5-S980 


15 


18 


16.42 


2 7 18.5 


14.264 


0.002 


46.88 


1.16 


4.525 


0.023 


-1.05 


0.01 



Table 3. This table lists the relevant data for the stars in the globular clusters NGC104, NGC3201, NGC4590, NGC5904 observed with 
VLT/FL AMES at LR, and that wc used in the analysis of the CaT-[Fe/H] calibration. For the CaT analysis, we select only those stars 
whose colour and magnitude are consistent with what is expected for RGB stars, and that have V — Vhb < 0, S/N> 10/A error in velocity 
< 5 km s -1 and velocities consistent with membership to the cluster (we assign membership to those stars within 3cr of the systemic 
velocity, see Table[TJ. The columns indicates: (1) the globular cluster name; (2) the star ID; (3), (4) star coordinates (right ascension in 
hours and declination in degrees); (5), (6) V magnitude and its error; (7), (8) heliocentric velocity and its error; (9), (10) summed CaT 
equivalent width (EW2 + EW3) and its e rror; (11 ) , (12) [Fe/H] value and its error, derived from the LR observations applying Eq. 1161 . 
The photometry and astrometry are from lStetsonl i2000l) . 
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0.07 
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2 39 4.31 


-34 25 18.9 
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0.004 


55.64 


1. 
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0. 
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-1 


.31 


0.07 


-0. 
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0. 
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-0.11 


0.06 
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-34 25 36.3 
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Table 4. This table lists the relevant data for the stars observed with VLT/FLAMES at both LR and HR resolution in the Sculptor and Fornax dSphs, and used in the analysis of the 
CaT-[Fe/H] calibration. The columns indicates: (1) the galaxy name; (2) the star ID; (3), (4) star coordinates (right ascension in hours and declination in degrees); (5), (6) V magnitude 
and its error; (7), (8) heliocentric velocity and its error; (9), (10) summed CaT equivalent width (EW2 + EW3) and its error; (11), (12) [Fe/H] value derived from the LR observations 
applying Eq. H16[l and its error; (13), (14) HR [Fe/H] va lue and error; (14), (15) [Ca/Fe] and its error, from the HR observations. The photometry and astrometry are from our ESO/WFI 
observations llTolstov et al. 2004; Battaglia et al. 200fJ). Those stars with [Ca/Fe]= —9.99 had too low signal-to-noise to allow a determination of the HR [Ca/H]. 



